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Abstract 

A measurement of the inclusive semileptonic branching ratios of b hadrons produced 
in Z decay is presented, using four million hadronic events collected by the ALEPH 
detector from 1991 to 1995. Electrons and muons are selected opposite to b-tagged 
hemispheres. Two different methods are explored to distinguish the contributions 
from direct b — * Xlv and cascade b —* c — * ~X.£v decays to the total lepton yield. One 
is based on the lepton transverse momentum spectrum, the other makes use of the 
correlation between the charge of the lepton and charge estimators built from tracks 
in the opposite hemisphere of the event. The latter method reduces the dependence 
on the modelling of semileptonic b decays. The results obtained by averaging the 
two techniques are 

BR(b->X&/) = 0.1070 ±0.0010 s t a t ±0.0023 syst ±0.0026 mode i , 
BR(b c Xlv) = 0.0818 ± 0.0015 stat ± 0.0022 syst t° Zl model • 
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1 Introduction 



Accurate knowledge of the direct inclusive semileptonic branching fraction of b hadrons, 
BR(b — > X£u) (l = eor/x), provides the opportunity to test and improve the theoretical 
models describing the dynamics of heavy hadrons. It is also an important input to 
a determination of \V c b\ Q, 0- Together with the cascade decay branching fraction 
BR(b —>■ c —>■ X£u), BR(b — > X£v) is an important input parameter for many heavy 
flavour analyses. 

Previous determinations of BR(b — > X£v) performed at the Z and the T(4S) have 
shown some disagreement, with that measured at the Z being higher while the opposite 
would be expected from the short b baryon lifetime ||. On the other hand, theoretical 
predictions have tended to be higher than experimental measurements although recent 
calculations, including higher order perturbative QCD corrections, give values in better 
agreement ||. ||. 

In this paper two new analyses are presented, based on the data collected with the 
ALEPH detector from 1991 to 1995. Two different methods are used to distinguish 
the contributions from the direct and cascade decays to the total lepton yield. One 
method has better statistical precision but a stronger dependence on the modelling 
of the semileptonic decays. The other is designed to have minimal decay modelling 
dependence. The efficiency of lepton identification is measured from data using several 
control samples. The description of the fragmentation of b quarks into b hadrons is based 
on the spectrum reconstructed with the ALEPH data and is therefore independent of 
modelling assumptions. 

2 The ALEPH detector 

A detailed description of the ALEPH detector can be found in |7j] and of its performance 
in ||. A brief overview will be given here, along with some basic information about lepton 
identification. 

Charged particles are detected in the central part of the apparatus which consists of 
a high resolution silicon vertex detector (VDET), a cylindrical drift chamber (ITC) and 
a large time projection chamber (TPC). 

The three tracking devices are immersed in a 1.5 T axial magnetic field provided 
by a superconducting solenoid; the particle momentum transverse to the beam axis is 
measured with a resolution cr (pr)/pr = 6 x 10~ A p T © 5 x 10~ 3 , with p^ measured in 
GeV/c. In the following, charged tracks are defined as "good" if a) they originate from 
within a cylinder, coaxial with the beam and centred around the nominal interaction 
point, with length 20 cm and radius 2 cm; b) their polar angle with respect to the beam 
satisfies the requirement |cos#| < 0.95; and c) at least four hits are reconstructed in the 
TPC. 

Hemispheres containing b quarks are identified with a lifetime b-tagging algorithm || 
which exploits the three-dimensional impact parameter resolution of the charged tracks. 
For tracks with two space points in the VDET (i.e. |cos#| < 0.7) the resolution can be 
parametrized as (25 + 95 /p) /mi, with p measured in GeV/c. 
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The TPC also provides up to 338 measurements of the specific ionization energy loss 
(dE/dx) which allows electrons to be separated from other charged particles by more than 
three standard deviations up to a momentum of 8 GeV/c. 

Electrons are identified by the characteristic longitudinal and transverse development 
of their associated showers in the 22 radiation length electromagnetic calorimeter (ECAL). 
The ECAL is segmented in 0.9° x 0.9° projective towers, read out in three longitudinal 
stacks; the energy resolution achieved for isolated electrons and photons is cr E /E = 
0.009 + 0. 18/ ^/E, with E measured in GeV. The dE/dx information provided by the TPC 
enhances the hadron rejection power, while non-prompt electrons originating from photon 
conversions in the material of the detector are rejected on the basis of their kinematic and 
geometric properties. 

Muons are identified by their penetration pattern in the hadronic calorimeter (HCAL); 
the additional three-dimensional coordinates measured in two layers of external muon 
chambers help in resolving the remaining possible ambiguities. 

The standard lepton identification technique is described in detail in [pTO . For the 



analyses presented in this paper, the selection is improved and optimized to reduce the 
systematic uncertainty on the results, as discussed in Section 01. 



3 Event samples 

The analysis is based on nearly four million hadronic Z decays selected using charged track 



information ]TT|. This data set was reprocessed during 1998 using improved reconstruction 
algorithms. The main benefits for this analysis are related to the enhanced particle 
identification capabilities, and are discussed in the description of the lepton selection 
(Section |). 

The statistics available for the simulation are larger than the data statistics by a 
factor of two; an additional sample of about five million Z — > bb simulated events is 
used for the estimate of the correlation between the b-tagging probability and the lepton 
selection efficiency (Section |6.14| ). Each event is divided into two hemispheres by the 
plane perpendicular to the thrust axis. Three samples of lepton candidates are selected 
as follows. 



Sample 3. A b-tagging variable £> tag , based on the large mass and lifetime of b hadrons, 
is built as in || (Fig. 0). This variable is defined using tracks contained in one hemisphere, 
but here the primary vertex is reconstructed using all tracks of the events, in contrast 
with || . The impact parameter in the simulation has been smeared in order to reproduce 
the resolution measured in the data, following the same procedure as in |J . The algorithm 
has good performance for events well contained in the vertex detector acceptance; events 
with | cos thrust | > 0.7 are rejected. 

A cut Stag > 2 is applied, selecting 345,555 hemispheres in the data, with a b purity 
of 97%. The cut is about 32% efficient on b hemispheres within the angular region 
considered. 

Lepton candidates (electron or muon) are searched for in hemispheres opposite to 
the selected ones. Events in which both hemispheres are b-tagged are used twice. The 



2 




Figure 1: Distribution of the b-tagging variable built by combining both lifetime and mass 
information of tracks in each hemisphere opposite to a lepton candidate. The error bars show 
the statistical error in the data only. 



candidates are ordered according to their transverse momentum p± , measured with respect 
to the jet axis as in |TUfl . When more than one lepton is found in a given hemisphere, the 
two leptons with highest p± and with opposite-charge are used for the analysis; if all have 
the same charge the one with the highest p± is taken. A total of 80,730 lepton candidates 
are finally selected. 



Sample T*. In each event one of the two hemispheres is randomly chosen, and a lepton 
candidate fulfilling the requirement p± > 1.25 GeV/c is searched for. This selects 148,001 
hemispheres in the data with an estimated b purity of 90% and a b efficiency of 12%. Since 
the p±_ cut suppresses cascade decays, the charge Q of the lepton is a good estimator of 
the charge of the quark at production time, leading to a probability of tagging the charge 
correctly of = 0.81 ± 0.01; this includes the dilution due to neutral B meson mixing. 



This value is measured from the data with the method described later in Section 5.4, and 



is in good agreement with the prediction of the simulation, P^(MC) = 0.817 ± 0.004 sta t- 
As this selection does not rely specifically on the vertex detector, no cut is made on 
| cos thrust |- The distribution of Q x p± is shown in Fig. |[ 

Leptons are searched for in hemispheres opposite to the tagged ones, as in the case of 
the B sample, yielding a total of 14,509 candidates. 
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Qxp ± (GeV/c) 

Figure 2: The charge estimator Q x p±, based on the charge of a high p± lepton: data are 
compared with Monte Carlo simulation. The contributions of hemispheres containing b and b 
quarks extracted from Monte Carlo are also shown. The cross-hatched area indicates the region 
rejected by the selection cut. 



The a priori random choice of the hemisphere used for charge tagging ensures that 
there is no double counting of lepton pairs, and allows the measurement of from the 
data (Section 



Sample Two hemisphere charge estimators are calculated: 



Qp 



s Ka 

where the sum runs over all the good charged tracks, as denned in Section ^ with 
momentum in excess of 200 MeV/c, q is the charge, p» the component of the momentum 
parallel to the thrust axis, s is the impact parameter significance, defined as in |§, n p = 0.5 
and k s = 0.3. Tracks with negative impact parameter are not included in the definition 
of Q s . 

The two charge estimators are combined using weights w, parametrised as a function 
of their magnitudes to achieve optimal performance on b hemispheres: 

Q H = w{\Q p \, \Q S \) Q p +[1- w(\Q p \, \Qs\)} Q s ■ 
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Qh 

Figure 3: Distribution of the hemisphere charge estimator Qh, m data and simulation. The 
contributions of hemispheres containing b and b quarks extracted from Monte Carlo are also 
shown. The cross-hatched area indicates the region rejected by the selection cut. 

Hemispheres are selected if they satisfy £> tag > 1.2, which enhances the b content of 
the sample. In addition a cut \Qh\ > 0.2 is applied to ensure a good probability that 
the sign of Qh is correlated with the charge of the b quark in the parent b hadron. 
The distribution of the Qh variable is shown in Fig. |3|, together with the contributions 
of hemispheres containing b and b quarks. Since the £> tag variable is used, events with 
Icos^thrustl > 0.7 are not considered. Hemispheres containing a lepton candidate with 
p± > 1.25 GeV/c are rejected in order to keep this sample statistically independent of 
sample V . 

The procedure selects 392,523 hemispheres in the data with an estimated b purity of 
87% and a b efficiency of 32%. The probability of correct b charge tagging measured from 
the data is pZ = 0.724 ±0.007, again in good agreement with the Monte Carlo prediction 
P b i(MC) =0.721 ±0.001 stat . 

The lepton yield in hemispheres opposite to the selected ones is studied as for the 
previous samples, giving a total of 61,012 candidates. 

4 Lepton identification 

A good control of the identification efficiency, as well as of the background in the selected 
sample, is crucial for this analysis. The identification of electrons and muons follows the 
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lines of [|T(J. A reduced dependence on the description of the b fragmentation is achieved 
by extending the acceptance to lower momentum leptons. 

The main change is, however, the use of a new estimator for the charged particle 
energy loss in the TPC. In the reprocessing of the LEP I data sample, information from 
the pulse height measured by the TPC pads has been used to build a new estimator 
similar to that from the wire measurements. The pad estimator is available for all tracks, 
while the wire estimator is calculated only for tracks that have a minimum number of 
isolated wire signals, which leads to an average inefficiency of about 15% in the hadronic 
environment. 

The wires estimator is defined as 

JW _ /JW\ . 

rw meas \ 11 

&? ' 

where (I w )i is the average energy loss expected for particle type i, of is the expected 
width, and /™ eas is the measured ionization. A similar equation defines the pad estimator 

If- 

A new energy loss estimator I is built from those two. This estimator is just the 
pad estimator for tracks that do not fulfil the requirement on the minimum number 
of isolated wire signals, otherwise it combines wire and pad information using the 
B.L.U.E. technique [0, taking into account the correlation of about 50% between the 
two measurements. 



4.1 Electrons 

For this analysis the momentum cut is lowered to p > 2 GeV/ c with respect to the selection 
described in ||10| . The requirement on the minimum number of isolated TPC wire signals 
is dropped, and a cut on the new energy loss estimator I e > —2 is applied. This removes 
the dependence of the identification efficiency upon the track isolation, and hence the 
electron p±. Compared to the previous selection, the new treatment of the energy loss 
information provides an increase in efficiency that goes from a few percent at high p±_ to 
about 30% at lowp^. The background increases marginally. 



4.2 Muons 

The momentum cut is lowered to p > 2.5 GeV/c. This still ensures that all muons reach 
the muon chambers but introduces a small dependence in the identification efficiency 
on the momentum. A dedicated study has been performed using 77 — > fi + fi~ events 
(Section 

A cut on the polar angle |cos#| < 0.7 is applied, which ensures that the muon is 
within the acceptance of both vertex detector layers. At least one VDET hit is required 
to be associated with the muon track, and the distance to the primary vertex in the plane 
orthogonal to the beam axis is required to satisfy \d Q \ < 2.5 mm. This substantially 
reduces the contamination from muons coming from kaon and pion decays. 

In addition a cut on the energy-loss estimator J M > —2 is applied; this further reduces 
the background from misidentified kaons by more than a factor of two. 
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5 The analysis method 



The lepton production rate in a sample of hemispheres with high b purity is dominated by 
the b — > X£V and b — > c — ■> X£z/ processes. In order to disentangle these two contributions, 
two different procedures are used, one based on the lepton kinematics, the other on the 
lepton charge. The two methods are described below. 

5.1 Transverse momentum analysis 

Sample £>, described in Section [| contains a large number of nearly unbiased b decays. 
The lepton rate in such a sample can be directly interpreted in terms of the sum 
of the direct and cascade inclusive semileptonic b decays, weighted by their selection 
efficiencies, once contributions from other physics sources and misidentified leptons have 
been corrected for. 

The study of the lepton rate as a function of the transverse momentum which 
discriminates the two components, allows the two branching ratios to be fitted 
simultaneously. 

Therefore a measurement of BR(b — > X£v) and BR(b —>■ c —>■ XJtv) can be obtained 
from a maximum likelihood fit to the number of leptons in each transverse momentum 
bin, as follows: 

where the product runs over the transverse momentum bins, rii is the number of leptons 
found in the data in each bin, \i{ is the expected number of leptons, which depends on the 
two branching ratios and contains the contributions from all the other sources of lepton 
candidates. 

Eq. p] can be rewritten as 

£ = x UHpi) , (2) 

counting s v 

p± spectrum 

where N is the total number of candidates observed in the data (N = J2i n i), A* i s the 
expected number (// = J^if^i), an d F{p±) is the binned function that gives the expected 
shape of the distribution of the candidates as a function of p± [Ti = fii/fi). The product 
runs over the lepton candidates. 

The expected number of candidates (/jl) in the selected hemispheres (A^emi) depends 
on the signal branching ratios multiplied by the average selection efficiencies e D ^ x , and 
on the rates of the other sources of lepton candidates: 

/i = iVhemi x [BR(b — > X£u)e\y^x£u + BR(b — > c — > X£z/)eb^ c -^x^i/ + other sources] . 

The expected p± shape depends on the shapes of the different components and their 
relative abundances: 
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where the shapes J ri (p±) are normalized to unity, the b — > X£u fraction is written as 

^ = BR(b -> X£u)e h ^ xeu 

BR(b — > X£z/)eb^x^ + BR(b — >■ c — > X£z/)eb-* c -+xfe/ + other sources ' 

and similar equations hold for all the components. 

The part of the likelihood labelled as "counting" contains the information on the total 
rate, and is therefore sensitive to the (weighted) sum of the two branching ratios. It is 
affected by uncertainties on the lepton identification efficiency and background, as well as 
on the rates of other processes yielding leptons. It has little dependence on the modelling 
of the p± spectrum. 

The part labelled as u p± spectrum" is sensitive to the relative contribution of the two 
signal sources, but almost insensitive to their overall normalization. It is heavily affected 
by uncertainties in the b decay modelling. 



5.2 Charge correlation analysis 

Another way of discriminating the b — > Xiu and b — > c — > X£u components is to exploit 
the different correlation of the lepton and the parent quark charges. The second part of 
the likelihood in Eq. |2| can be replaced with a term containing the fraction of leptons that 
have the same or the opposite charge relative to a charge estimator built using tracks in 
the opposite hemisphere (e.g. jet charge). 

However, such a method would have poor statistical power if applied to sample B. 
The charge tag samples V and J are selected for this reason, somewhat relaxing the 
requirement on the b purity in favour of higher statistics and better charge tagging. 

A likelihood function is constructed based on the counting part from sample B and 
the charge spectra of samples V and J as follows: 

C= "-^f x (1-F v f* x J^'jl-Fj)^ , (3) 

V 7" TL\ charge (V) charge (J) 

counting (a) 

where Tj> is the expected fraction of lepton candidates with charge opposite to the charge 
tag of the other hemisphere in sample V, Nq v and N$ v are the numbers of candidates 
with opposite and same charge found in the data. The same holds for sample J . 

The expected fractions T-p and Tj are sensitive to the relative contributions of the 
b — > X£u and b — > c — > Xiu components and depend on the rate of correct tagging for 
the opposite hemisphere charge estimator, as well as on the background components. 



5.3 Flavour composition of the selected samples 

The flavour composition for each of the three samples is estimated as follows. The fraction 
-Phemi of single b-tagged hemispheres is measured from the data. The efficiency e c for 
tagging charm events and the average efficiency e u d s for tagging light quark events are 
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determined using simulated events, and the sample composition is calculated as 

„b _ -I _ Rc^c + (1 ~ -Rb ~ Pc) e uds 

J hemi p ) 

-fnemi 

R c e c 

J hemi p ' 

-fnemi 

(•uds (1 R& -Rc) ^uds 

J hemi p ' 

-fnemi 

where Pb and R c [|13j are the ratios of the bb and cc partial widths of the Z to the total 
hadronic width, taken from experimental measurements. 

5.4 Charge tagging in samples "P and 3 

The terms T-p and Tj in Eq. |3] are written in terms of the probabilities that the charge of 
the parent quark in the hemisphere opposite to the lepton candidate is correctly tagged. 
These probabilities for b or b quarks P bb , discussed in Section [|, are measured from the 
data by means of a double-tag method. The description below applies to both samples V 
and J . 

The charge-tagging selection cut is applied to both hemispheres, and the events 
selected are split into three classes: those in which the two hemispheres have opposite 
charge, and those in which have the same charge, positive or negative. The three fractions 
P oc , F ++ and F can be written as 

poc e b poc _i_ fc poc , i-uds poc (a \ 

r ~ Jevt r b ' Jevt r c ' J evt r uds > 

where f^ vt , f£ vt and are the contributions of b, c and light flavour events to this 
sample, with similar equations for F ++ and F . The different flavour contributions are 
measured with the same procedure used to estimate the hemisphere sample composition 
(Section |573|) : 



„b RcCc + (1 — Rb — P C )Cuds 

J evt rn 

-Tevt 

, c = RcCc 
J evt T-i ' 

-^evt 

„ uds = (1 ~ -Rb - -Rc) Cuds 
J evt t-i ' 



-c'evt 

where P ev t is the fraction of events with both hemispheres tagged in the data, and Cc and 
Cuds are the charm and light quark event efficiencies determined from the simulation. 

The fractions F° c (P c ++ , F c ) and F°^ s (F^, F~^) of events with opposite (same) 
hemisphere charges in charm and light quark events are taken from the simulation, and 
Eq. | is solved for P b oc (P++, F b "~). 

In order to extract P bb from P b c , several effects must be considered. 

• Because of the interaction of particles with the detector material, the probability for 
correct charge tagging is different for b and b quarks. The difference can be defined 
as P b DT = P b ^ T (l - e), P b DT = P b I b T (l + e), where P b ^ T is the average probability of 
correct charge tagging in the double-tagged sample selected. 



9 



• The probabilities of correct charge tagging in the two hemispheres are not 
independent. The correlation p is defined so that a hemisphere opposite to one in 
which the charge was correctly tagged has the probability of correct charge tagging 
enhanced by a factor (1 + p). 

• The sample used for the branching ratio measurement includes hemispheres for 
which the opposite one did not fulfil the selection cuts. The selection in the 
opposite hemisphere has some correlation with the charge tagging probability of 
the hemisphere considered, therefore the parameter which is measured with the 
double tag method requires a further correction, a, as explained below. 

With all these effects taken into account, -PJr T is related to P b c via 

K = (^ T ) 2 + (l - ^ T ) 2 + 2(P b f) 2 (P - e 2 - e 2 p ) , (6) 
where e can be written as 

PI++ - PI- 
CT) 



2^ T 



Eq. |] and ^ can be solved for P^ T and e taking for p the estimate provided by the 
simulation. Then P bb is calculated as 



DT 



P bb = (1 - a)P h „ 

in terms of a bias correction a also taken from simulated events. Both a and p are found 
to be smaller than 1%. 



6 Sources of systematic errors 

In this section the sources of possible systematic effects are discussed. The estimated 
uncertainties are summarized in Tables [I] and 0. 

6.1 Z partial widths to bb and cc 

The values of P b and R c are used in the derivation of the sample compositions of the three 
hemisphere samples, and in the calculation of P^ and P^. The most recent averages are 
taken 0, P b = 0.21653 ±0.00069 and P c = 0.1709 ±0.0034, and the experimental errors 



are considered as sources of systematic uncertainty. 



6.2 Heavy quarks from gluon splitting 

Charm and bottom quark pairs may be produced from the splitting of hard gluons. 
The heavy flavour hadrons resulting from this process have a significantly softer energy 
spectrum and thus give rise to a source of prompt leptons with kinematic properties 
substantially different from those produced by heavy hadrons from direct Z decay. In 
addition, leptons originating from gluons splitting to heavy quarks have a random charge 
correlation with the charge estimators defined in the opposite hemisphere. 
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A[BR(b 

L \ 




A[BR(b -> 


c -> X£u)\ 


Source 


P± 


Charge 


Px 


Charge 




negl. 


negl. 


negl. 


negl. 


Rc 


± 0.005 


T 0.007 


T 0.002 


± 0.017 


N(g -> bb) 


T 0.002 


T 0.002 


T 0.002 


=F 0.001 


N(g -» cc) 


T 0.001 


T 0.006 


T 0.014 


=F 0.006 


electron ID efficiency 


T 0.063 


T 0.081 


T 0.087 


T 0.056 


7 conversions 


± 0.003 


T 0.006 


=F 0.022 


T 0.008 


electron bkg 


± 0.004 


T 0.007 


=F 0.026 


=F 0.009 


muon ID efficiency 


± 0.065 


± 0.063 


± 0.039 


± 0.039 


muon bkg 


± 0.002 


T 0.013 


T 0.037 


T 0.015 


BR(b -> c X£i/)|^ 


± 0.004 


± 0.022 


± 0.002 


T 0.026 


BR(b -» J/^ (^') -» #) 


negl. 


negl. 


negl. 


negl. 


BR(b -» r -» £) 


T 0.017 


T 0.043 


=F 0.053 


T 0.011 


BR(b -> W -»■ c -» ^) 


± 0.010 


T 0.223 


=F 0.407 


T 0.039 


BR(c X&/) 


negl. 


=F 0.016 


T 0.009 


± 0.016 


BR(b -> X u £) 


T 0.032 


=F 0.022 


± 0.013 


T 0.004 


b fragmentation 


T 0.074 


=F 0.089 


T 0.120 


T 0.101 


c fragmentation 


± 0.001 


± 0.005 


negl. 


T 0.005 


e c sample B 


± 0.027 


± 0.015 


T 0.009 


=F 0.010 


e u ds sample S 


± 0.015 


± 0.016 


± 0.012 


± 0.011 


e c sample J7 


_ 


T 0.018 


_ 


± 0.029 


e uds sample J 


_ 


negl. 


_ 


negl. 


e c sample P 


_ 


T 0.012 


_ 


± 0.019 


e uc j s sample V 


_ 


negl. 


_ 


negl. 


c charge tag rate 


- 


± 0.036 


- 


T 0.057 


b charge tag rate 


- 


± 0.069 


- 


T 0.109 


Mixing in b — * Xiv 




± 0.035 




T 0.055 


Mixing in b — > c — > Xlv 




T 0.055 




± 0.087 


bkg charge correlation 




± 0.027 




T 0.043 


b tag - lept correlation 


± 0.006 


T 0.007 


T 0.025 


T 0.005 


Total 


± 0.128 


± 0.290 


± 0.443 


± 0.212 



Table 1: Estimated contributions to the systematic uncertainties on BR(b — > X£zv) and 
BR(b — ► c — > X£i/)). Results for both transverse momentum and charge correlation analyses 
are given. Uncertainties related to the modelling of semileptonic decays are shown separately in 
Table ||. All values are given in units of 10 -2 . 



The latest world average values are used for the number of gluons splitting to heavy 
quarks per hadronic Z decay p4| , 



iV(g-+bb) = 0.00254 ± 0.00051 , 
iV(g -> cc) = 0.0296 ± 0.038 , 

and the experimental errors are used to estimate the associated uncertainty. 
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A[BR(b 


yj 


A 


BR(b -> 


c — > 


X£v)} 


Source 


P± 


Charge 


P± 


Charge 


Y)lv 


± 0.012 


± 0.022 


± 


0.029 


± 


0.013 




T 0.008 


± 0.077 


± 


0.035 


± 


0.010 


Inclusive D^Xlu 


± 0.254 


± 0.086 




0.265 


=F 


0.004 


\J\IV 


=F 0.068 


=F 0.028 


± 


0.056 


=F 


0.007 




=F 0.047 


T 0.018 


± 


0.040 


=F 


0.006 


Broad states 


± 0.327 


± 0.169 




0.191 


± 


0.052 


Bg fraction 


± 0.049 


± 0.060 


± 


0.063 




0.048 


Baryon fraction 


T 0.025 


=F 0.012 


± 


0.059 


± 


0.044 


b — ► X£i/ modelling 


± 0.426 


± 0.202 


± 


0.348 


± 


0.085 


c — > ~X.£v modelling 


-0.087 
+0.072 


-0.038 
+0.021 




-0.037 
-0.020 


-i 


0.117 
-0.063 


b — > D modelling 


-0.072 
+0.060 


-0.002 
+0.001 




-0.055 
-0.049 


-t 


-0.058 
-0.050 


Total modelling 


+0.436 
-0.441 


+0.204 
-0.206 


4 


-0.353 
-0.355 


-i 


-0.117 
0.156 



Table 2: Estimated contributions to the systematic uncertainties on BR(b — > X£^) and 
BR(b — > c — > X£z^)) related to the modelling of semileptonic decays. The uncertainty in the 
modelling of direct b — ► Xlv decays is estimated by adding in quadrature the uncertainties 
from the rate of the different charmed mesons, inflated by 25%, and the uncertainties from the 
production rates of Bg mesons and b baryons, as discussed in Section 5.15. Results for both 



transverse momentum and charge correlation analyses are shown. All values are given in units 
of 10" 2 . 



6.3 Electron identification efficiency and background 

The electron identification efficiency is measured from data using photon conversions in 
the detector material. Correction factors are derived, with respect to the Monte Carlo, 
for the dependence on momentum, transverse momentum and polar angle. These factors 
typically differ from unity by less than 1%. The associated systematic uncertainty is 
estimated by removing the corrections. 

The background from hadron misidentification is estimated by removing the cut on 
the energy loss, and studying the shape of the I estimator. This study is performed on 
data and Monte Carlo and no significant deviation is observed. An uncertainty of 20% 
on the backgroud level is assigned from the statistical precision of the method. 

The background from unidentified photon conversions is estimated by studying the 
shape of the variable 

p 7 =qd pT , 

where q is the charge of the electron, do is the distance of closest approach to the primary 
vertex in the xy plane signed according to the track angular momentum around the 
origin, and pp is the component of the track momentum transverse to the beam axis. The 
variable p 7 is expected to be zero for prompt electrons and proportional to the square of 
the materialization radius for electrons coming from photon conversions. The study of 
the positive tail of the distribution yields a correction factor of 1.05 to be applied to the 



12 



simulation, with a statistical error of 0.02. The correction factor is removed to estimate 
the systematic uncertainty. 

6.4 Muon identification efficiency and background 

The identification efficiency for high energy muons is measured from the data using Z 
decays to muon pairs, as a function of polar and azimuthal angle. Simulated events are 
reweighted to reproduce the measured efficiencies. Correction factors are typically a few 
per mil. 

Simulated events show some dependence of the identification efficiency upon the muon 
momentum for momenta around 3 GeV/c. This effect is also checked on real data using 
77 — > events. Additional correction factors, of order 1 to 2% are derived for muons 

with momenta between 2.5 GeV/c and 4 GeV/c. The shift observed when the corrections 
factors are removed is used as an estimate of the systematic uncertainty. 

The main background for muon candidates is due to misidentified pions as well as pions 
decaying before entering the calorimeters. The corresponding contributions for kaons are 
substantially reduced by the cut on the measured energy loss, and are estimated to be a 
factor of four smaller. 

In order to check the background rate from data, Kg — > 7r + 7r~ decays are selected in 
hadronic events to yield a 99% pure sample of pions. The muon identification procedure is 
applied to these tracks, and the selection efficiency is compared between data and Monte 
Carlo. Agreement is found within the statistical precision of the test, which is 5%. The 
test is repeated applying different £> tag cuts in the hemisphere opposite to the Kg candidate 
in order to check for a possible dependence on the flavour. No trend is observed. 

The uncertainty of 5% estimated from the check with Kg — > 7r + 7r~ is enlarged to 
10% for the assignment of a systematic error to the muon background. This allows for 
additional uncertainties from the smaller kaon component, as well as possible differences 
between data and simulation in the production rates and kinematic properties of pions 
and kaons in b events. 

6.5 Specific b c — > X£v correction 

The b — > c — > X£v rate is in principle different in transitions where the W from the 
b hadron decays leptonically, BR(b — > c — > X£u)\ e , or hadronically, BR(b — > c — > X£v)\ h . 
The fit could yield a biased result for the average BR(b — > c — > X£v) if the acceptance 
were different for the two cases. This effect is investigated by changing in the simulation 
the relative population of the two species and recalculating all efficiencies and spectra. 
The BR(b — > c — > X£u)\ e is increased to a conservative 20% and the BR(b — > c — > X£v)\ h 
is decreased accordingly. The shift observed in the fitted values is taken as an estimate 
of the associated systematic error. 
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6.6 Other sources of prompt leptons 

The rates of leptons coming from J/0 and from intermediate r decays used for this analysis 



are 



BR(b -> J/^ (VO -> ^) = 0.00072 ± 0.00006 , 
BR(b -> r -> £) = 0.00419 ± 0.00055 . 

Leptons produced from cascade b decays where the intermediate charm is produced 
from a W — > cs transition, denoted b — > W — > c — > £, are also a background to the 
analysis. They affect most directly the result for BR(b — > c — > X£u) in the transverse 
momentum analysis, as they have kinematic properties similar to b — > c — > X£z/ decays. 
On the other hand in the charge correlation analysis only the value of BR(b — > X£z/) 
depends on the rate of these transitions since the correlation between the charge of the 
lepton and the charge of the parent quark is the same as in b — > X£z/ direct decays. The 



value used is [|14[] 

BR(b -> W -> c -> £) = 0.0162 ± 0.0044 . 

The semileptonic decays of charmed hadrons in the residual background of cc events 
contribute to the total observed lepton yield. The LEP average value |T3 is taken: 



BR(c i) = 0.0984 ± 0.0032 . 

All the above sources of leptons are subtracted from the total lepton yield. 

The rate of lepton production in charmless semileptonic b decays, b — > X u £, mainly 
affects the high end of the p± spectrum. A variation of this contribution relative to the 
total BR(b — > X£v) is included in the systematic error calculation. The value used in this 
analysis is JOJ 

BR(b -> XJ) = 0.00167 ± 0.00055 . 



6.7 Fragmentation of b quarks 

The scaled energy spectrum of b hadrons in the simulation is modified in order to 
reproduce the spectrum reconstructed in the model-independent analysis of 0. 

The statistical and systematic uncertainties on the population of each energy bin are 
propagated to the measured branching ratios, taking into account bin-to-bin correlations. 
The systematic errors on the energy spectrum due to the uncertainty on the charmed 
meson species produced in B meson decays are not considered here, as they are correlated 
with the uncertainty on the modelling of semileptonic decays (Section |6.15| ). 

The correction factors derived from the comparison of the measured and simulated 
energy spectra of B mesons are also applied to B° mesons and b baryons. 



6.8 Fragmentation of c quarks 



Charm fragmentation is simulated using the Peterson et al. model |L6[. The parameter 
e c , which controls the shape of the function, is adjusted to reproduce the measured value 
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of the mean scaled energy of weakly-decaying charmed hadrons, (x c ) = 0.484 ± 0.008 fl4 |, 
and then varied by its uncertainty. 



6.9 Charm and light quark background 



The calculation of the sample compositions described in Section |5.3| relies on the 
simulation for estimates of the charm and light quark hemisphere selection efficiencies. 
For samples B and J the estimate of the uncertainty on the background efficiencies follows 
the lines of |J, where all the relevant physics parameters in the Monte Carlo have been 
reweighted. The values used for this analysis are 

e f = 0.00939 ± 0.00094 e^ ds = 0.00060 ± 0.00015 , 
ef = 0.0439 ±0.0022 ef ds = 0.0054 ±0.0008 . 



For sample V the dominant sources of uncertainty are the charm semileptonic 
branching fraction and decay modelling, for charm hemispheres, and lepton background 
for light quark hemispheres. The corresponding values and errors are 

e v c = 0.0108 ±0.0011 e^ ds = 0.0016 ± 0.0003 . 

6.10 Charge tagging in charm and light quark hemispheres 

The fractions of opposite charge charm and light quark hemispheres entering Eq. ^ are 
taken from the simulation. The difference between F ++ and F for these flavours is 
neglected in the analysis, therefore they are determined once F oc is obtained. 

In sample V, F° c is determined by the relative amounts of true and fake leptons 
selected, and by the degree of correlation between the charge of fake leptons and the charge 
of the c quark. The rates of lepton production from c hadron semileptonic decays and from 



hadron misidentification are varied within their estimated uncertainties (Sections |6.3| , pA 
and |6.6| ). The charge correlation for fake leptons is set to zero and half of the difference 
in F° c is taken as an estimate of the related uncertainties. Adding the uncertainties from 
the three sources in quadrature yields F c oc = 0.87 ± 0.04 in sample V, corresponding to 
= 0.93 ± 0.02. 

For sample J an estimate of the uncertainty on F c oc is obtained by comparing the 
value with Fg°. For charm hemispheres, the correlation between the charge of the quark 
and the sign of the estimator is stronger than in the case of b hemispheres. This is 
valid separately for the two estimators Q p and Q s . For the first, the reason is the higher 
charm quark charge. In the case of Q s , the effect is due to the b-tagging cut which 
selects charm hemispheres with particularly long lifetime, and to the smaller charged 
particle multiplicity of c hadron decays compared to b hadron decays. The systematic 
uncertainty is evaluated by setting the charge correlation for charm hemispheres equal to 
that for b hemispheres and taking half of the difference as the systematic uncertainty; 
this procedure yields F° c = 0.64 ± 0.02 in sample J7, corresponding to P£ = 0.77 ± 0.02. 
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In light quark events, F°^ s is still somewhat larger than 0.5 due to the correlation of 
the estimator with the parent quark charge, for both samples. However the systematic 
uncertainty obtained by setting it equal to 0.5 is negligible. 



6.11 Charge tagging in b hemispheres 

Besides the sources already considered, the uncertainty on P bb also depends on the 
uncertainties on P oc , P ++ , F , F evt , ( c , £ u d s , p and a appearing in Eq. f|-§. 

The statistical errors on F oc , F ++ , F and P ev t, which are measured from the 
data, are propagated to the results of the fit and included in the statistical errors on 
the branching ratios. 

For sample J the values of p and a measured in the simulated events are 0.0096 and 
0.0063, respectively. For sample V, a is found to be 0.0062, while there is no significant 
correlation between the two hemispheres. The associated systematic uncertainty is 
estimated by setting p and a to zero and taking half of the observed shift. 

The uncertainties from ( c and Cuds are negligible. 



6.12 Neutral B meson mixing 

The mixing of neutral B mesons contributes to the degradation of the correlation between 
the charge of the lepton and the charge of the parent b quark produced in the Z decay. 
The LEP average \ = 0.1194 ± 0.0043 jn§ is used as input in the likelihood. This value 



is interpreted as the average mixing rate for the b hadron mixture from b — > X£z/ decays. 

The relative population of B° and B + mesons is not equal in b — > c — > X£z/ decays, 
due to the different semileptonic branching ratios of D + and D° mesons, leading to 
an effectively higher value of the average mixing parameter = % (1 + <5). The value 
5 = 0.13 is estimated from the simulation and is varied by 50% to estimate the systematic 
uncertainty. 

As a consequence the expected fraction of events with hemispheres of opposite sign 
for the b — > X£u and b — > c — > X£z/ components is 



X£u: T = P b b(l-X) + (1- Ab)x 
X£u: T = P bb x / +(l-P bB )(l-x 1 



-,/\ 



6.13 Charge correlation for lepton background 

Leptons coming from kaon and pion decays in flight as well as misidentified kaons and 
pions retain some information about the charge of the primary quark, both in charm 
and in b events. This effect must be taken into account when evaluating the opposite- 
charge and same-charge fractions in the charge correlation analysis. The rate at which 
the information about the quark charge is retained is measured in the Monte Carlo. The 
systematic uncertainty is evaluated by taking half of the difference between the Monte 
Carlo value and no charge correlation. 
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6.14 Hemisphere-hemisphere correlation 

The b-tagging efficiency is affected by the presence of a lepton in the hemisphere opposite 
to that used for the flavour tagging, in a way that depends on both the momentum and 
the transverse momentum of the lepton. In particular, hemispheres opposite to a high 
momentum lepton are more likely to satisfy the b-tagging cut since such leptons are likely 
to be produced in events with little gluon emission, where both b hadrons have high 
momentum. The b-tagging efficiency for those hemispheres is then higher than expected 
in an unbiased sample of hemispheres. 

The ratio R = e b /(eb x e{) is determined from simulated events, where eb and e\ 
are respectively the b-tagging efficiency and the probability of finding a lepton in a 
b hemisphere while e b is the probability that in the same event one hemisphere is fa- 
tagged and in the opposite hemisphere one lepton is selected. The dependence of R 
upon the lepton momentum is irrelevant for the analysis, while the overall value affects 
the part of the likelihood labelled as "counting" in Eq. § and |], and the dependence on 
the lepton transverse momentum affects the part labelled as u p± spectrum" in Eq. |2|. 
The quantity R is therefore determined as a function of the lepton p±, for the different 
sources of lepton candidates, and the corresponding correction factors are applied in 
the analysis. Deviations from unity are typically smaller than 1%. The corresponding 
systematic uncertainty is estimated by removing the corrections. 



6.15 Modelling of direct semileptonic b decays 

In the past, in all measurements of the inclusive b — * Xiu rate the decay kinematics 
were modelled using the ACCMM spectrum, and the associated systematic uncertainty 
was estimated from the shift observed using the ISGW and ISGW** spectra, as proposed 

0- 



in 



In this work a different approach is used, based on the available measurements of the 
B meson decay rates into the different charm hadron final states. B° and B + mesons 
decay semileptonically into D, D* and D** mesons, as well as non-resonant D^7T final 
states. Leptons coming from each of these components have a different energy spectrum 
so that the shape of the inclusive B — > Xiv spectrum depends on the branching fractions 
of B° and B + into the various charmed species. 

The fractions of D, D*, Di, and D£ in the simulation are reweighted using the latest 
measured values |15|, |18| reported in Table ||. The broad D** states are assumed to be 



equal to the sum of the narrow D x and D£ states, and the non-resonant D^*'tt decays 
account for the rate needed in order to add up to the measured inclusive D^X branching 
ratio. 

The systematic uncertainty is evaluated by varying the measured branching fractions 
within their estimated errors. Additionally the rate of the broad states is set to zero 
and compensated entirely with the non- resonant D^ir states. In each case, the B meson 
energy spectrum measured in |J is used, thus taking into account the correlation with 
the present analysis. This approach renders the analysis presented independent of the 
modelling of the b fragmentation in the simulation. 

The energy spectra obtained with this procedure for the mixture of electrons and 



17 



Process 


BR (%) 


B -> Dlv 
B -> D*£i/ 
B -> DWXfe 
with B -> Di^i/ 
with B -> D*£z/ 


1.95 ± 0.27 
5.05 ±0.25 
2.7 ±0.7 
0.63 ±0.11 
0.23 ±0.09 



Table 3: Branching ratios for semileptonic decays of B mesons with different charmed mesons 
in the final state. 



muons used in the analysis are compared in Fig. |] with the inclusive spectra given in 
In Fig. [5] the effect of the corrections applied to estimate the systematic uncertainty is 
shown. The spectra resulting from this semi-exclusive treatment always lie between the 
softest (ISGW**) and hardest (ISGW) spectra. 

The procedure described applies to B° and B + decays, which represent about 80% of 
the sample. Uncertainties in the production rates of B° mesons and b baryons can affect 
the analysis to the extent that their decay kinematics differ from those of B° and B + 
mesons. The B° and b baryon production rates are reweighted to the latest experimental 
results, f B o = 0.094 ± 0.022 and f b -bar y on = 0.101 ± 0.018 PH. Uncertainties in the 



modelling of B° and b baryon decays are accounted for by enlarging the uncertainties from 
the branching ratios of Table |]by 25%. Finally, the observed shifts on BR(b — > X£v) and 
BR(b — > c — > XEis) are symmetrized and added in quadrature. 



6.16 Modelling of c — > X£z/ and b — > c — ► X£v decays 



The c — > X£z/ decay spectrum is obtained by means of a combined fit to DELCO [19| and 
MARK III PDA data, performed using the ACCMM model and varied as described in [17 . 

In order to describe b — > c — > X£v decays, the model proposed for c — > Xiu is combined 
with the measured b — > D spectrum from CLEO [21] following the procedure of [17]. 



7 Results 

The results obtained with the transverse momentum analysis are 

BR(b -> yjLv) = 0.1107 ± 0.0007 stat ± 0.0013 syst ± 0.0044 model , 
BR(b -> c -f X£v) = 0.0752 ± 0.0010 stat ± 0.0044 syst ± 0.0035 mode i , 

with a statistical correlation of —0.45. The fitted spectrum is compared to the p± 
distribution observed in the data in Fig. [| where the contributions of the two signal 
processes are also shown. 

The charge correlation analysis gives 



BR(b-fX^) = 0.1057 ±0.0011 stat ±0.0029 syst ±0.0020 modc i 

syst -0.0016 model , 



BR(b -> c -> YJLv) = 0.0830 ± 0.0016 stat ± 0.0021 -"■ " 1 ' 2 
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Figure 4: Lepton energy spectrum in the b hadron rest frame. The histograms show the 
distributions obtained after reweighting with the ACCMM, ISGW and ISGW** models. The 
dots show the spectrum after correcting the Monte Carlo by following the procedure described 
in Section 3.15 using the central value of the branching ratios listed in Table 13. 



with a statistical correlation of —0.76. Since this method does not use the information 
from the lepton p±, the p± spectrum can be extracted from the data by repeating 
the charge correlation analysis in each p± bin. The values of BR(b — > X.£u) and 
BR(b — > c — > X£v) measured in each p± bin are then used to reweight the Monte Carlo 
spectrum previously corrected as described in Section |6.15| . In Fig. [7| the b — > X£v and 
b — > c — > X£u spectra shapes from the data are compared with the Monte Carlo. The 
discrepancy observed reflects the difference between the results obtained with the two 
methods. 

The differences can be quantified taking into account the correlation between 
systematic errors, as A b _o«K = 0.0049 ± 0.0036 (total error, of which 0.0025 is due 
to modelling) and A b _> c _o«K = 0.0077 ± 0.0056 (total error, of which 0.0039 is due to 
modelling). The two results are therefore consistent within 1.4 times the total uncorrelated 
error. 

The statistical errors are correlated and smaller than the systematic errors. In the 
calculation of the average, maximal statistical correlation is assumed between the two 
measurements, i.e. 81% for BR(b — > X.£u) and 85% for BR(b — > c — ► X£v), which 
leads to a slight overestimate of the statistical errors. The weights used to calculate 
the average are W b -+xfe(p_i_)= 0.25, W b -oc&, (charge) = 0.75, W b ^ c _ X fe(p±)= 0.15, 
W b ^ c _»xfe (charge) = 0.85. They are determined by minimizing the total error (B.L.U.E. 
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Figure 5: The plots (a) and (b) show the ratio between the lepton spectra given by the ISGW 
and ISGW** models and the spectrum obtained by correcting the Monte Carlo using the central 
value of the branching ratios listed in Table [||. Histograms (c)-(h) show the effect on the 
spectrum due to the variation by one standard deviation of each component with respect to the 
central value. 



technique [fL^Q. The resulting average values are 

BR(b->X£z/) = 0.1070 ±0.0010 stat ±0.0023 syst ±0.0026 mode i 
BR(b -> c -> X£u) = 0.0818 ± 0.0015 stat ± 0.0022 syst +°;°°}° model . 



7.1 Consistency checks 

The analysis was performed seperately for each year of data taking, to check for possible 
detector-related effects. The results are consistent throughout the whole LEP I period. 
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Figure 6: Comparison between the lepton pj_ spectrum measured in the data and the result of 
the fit. The contribution from b — > X£v as well as the sum of b — > X£z^ and b — > c — ► X&/ are 
shown. The error bars show statistical errors on the data. 



The analysis was also performed using electron and muon candidates separately. The two 
species of leptons have different uncertainties in the selection efficiency and in the rate of 
background. The tighter selection cut applied on the muon momentum is reflected in a 
different dependence on modelling. The results are consistent within errors as shown by 
their differences quoted below. 

Transverse momentum analysis: 

K'Sxeu = -0.0029 ±0.0013 stat ± 0.0018 syst ± 0.0002 modcl , 
= 0.0017 ± 0.0020 stat ± 0.0021 syst ± 0.0011 modcl . 

Charge correlation analysis: 

K'Sxe, = -0.0015 ±0.0021 stat ± 0.0021 syst ± 0.0004 modcl , 
Ablfc^ = -0.0009 ±0.0031 stat ±0.0017 syst ± 0.0005 modcl . 

The cuts applied on the b-tagging variable, on the jet-charge value and on the p± were 
also varied. The results, in both methods, were found to be stable. 
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Figure 7: Comparison between the lepton p±_ spectra measured in the data with the charge 
correlation analysis and that expected from simulation, (a) b — * lL£v\ data are compared 
with the simulation corrected as described in Section 6,15| , ISGW and ISGW** models, (b) 
b — > c — > Xlv. All the spectra are normalized to unit area. 



8 Conclusions 



Two different strategies have been adopted to distinguish the b — * X£u and the 
b — > c — > X£z/ decays in order to measure the inclusive semileptonic branching fractions 
of b hadrons, BR(b — > X£v) and BR(b — > c — > X£u). One measurement, based 
on the analysis of the lepton transverse momentum, is significantly affected by the 
uncertainty from the modelling of the b — > X£v p± spectrum. The other is based 
on the charge correlation between the lepton and a charge estimator in the opposite 
hemisphere. This has smaller dependence on the modelling of the decay kinematics, but 
the determination of the b — > X£v branching fraction suffers from the uncertainty on the 
rate ofb— > W —»■ c ^ £ decays. In both cases the modelling of the inclusive b — > Xlv 
decay kinematics is based on the measured rates of the different c hadrons in the final 
state. 

The values of BR(b — > X£u) and BR(b —>■ c — > X.£v) measured with the two methods 
are consistent within errors and they have been averaged. The resulting values are 



BR(b 



BR(b 



X£u) 
X£v) 



0.1070 ±0.0010 stat ±0.0023 syst ±0.0026 modcl 



0.0818 ±0.0015 stat ± 0.0022 syst ^;™J model , 



in agreement with the latest measurements performed by other LEP experiments [^, |23 



24], [25], and with the latest values from experiments running at the T(4S) resonance |T8 
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